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the electrical characteristics of the system. However, while there is a wide experience in modeling the 
behavior of traditional photovoltaic systems, not every model for flat-plate solar cells or modules is 
directly applicable to high concentration photovoltaic cells or modules because of the special features of 
these devices (use of multijunction cells, use of optics for high concentration, etc.). So, in recent years, the 


Keywords: scientific community has devoted considerable efforts in developing models that reproduce the electrical 
Concentration photovoltaics behavior of high concentration cells and modules. These models allow calculating the main electrical 
Cell model parameters of the device from its operating conditions (irradiance, cell temperature, spectral distribution 


Module model of the radiation, etc.). In this paper, a comprehensive review of existing models for the electrical 


characterization of high concentration photovoltaic cells and modules is presented with the aim of 
helping the photovoltaic professionals and researchers in the design, monitoring and energy 
prediction tasks. 
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1. Introduction 


High Concentration PhotoVoltaic (HCPV) cells and modules oper- 
ate under concentration ratios between 300 and 2000 suns. The 
concentration mechanism is done by lenses or mirrors that either 
reflect or refract the sunlight on the solar cell surface. Replacing the 
expensive solar cell materials by cheaper optics offers the opportu- 
nity of reducing the costs of the modules. In addition, the use of 
multijunction solar cells implies an increase of the system efficiency, 
what also can reduce the associated costs [1]. 

Three main elements are required in order to generate elec- 
tricity: the multijunction solar cell, the optics and the solar tracker. 
These elements can be combined in a lot of different configura- 
tions [2,3]. A typical commercial HCPV system is composed of 
concentration modules mounted on two-axis solar trackers, where 
each module contains several multijunction solar cells and each 
cell has its own optics, typically a point-focus Fresnel lens. Triple- 
junction GalnP/GalnAs/Ge cells are usually utilized [4], while other 
material combinations are now under research [5-8]. Other 
auxiliary elements are used, such as heat sinks, Secondary Optical 
Elements (SOEs), bypass diodes and protection elements. 
Nowadays, these configurations dominate the HCPV market and 
are very promising solutions for the future of photovoltaics [9]. 

Modeling the electrical behavior of HCPV cells and modules is 
an important task for the design, monitoring and energy predic- 
tion of these kinds of systems. In this paper, a review of the state- 
of-the-art of models that reproduce the operation of HCPV cells 
and modules is presented. These models calculate the I-V char- 
acteristic curve of the device from its operating conditions. Not 
every model allows calculating the whole I-V curve; some of them 
only calculate some relevant electrical parameters, such as the 
maximum power (Pypp), the short-circuit current (Isc), the open- 
circuit voltage (Voc), the fill factor (FF), etc. 

The presented models differ in the parameters required to 
characterize the device, in the variables used as inputs and, finally, 
in the equations and procedures used to obtain the electrical 
characteristics. The models present different levels of complexity 
and accuracy and require different equipments to be applied. This 
review remarks the characteristics of each one of these models in 
order to help in the choice of the most suitable model for each 
specific application. 

The paper is organized as follows: in the second section, a brief 
overview of the technologies under study is presented; in the third 
section, the existing models for multijunction solar cells are 
summarized; the fourth section describes the main contributions 
in the modeling of HCPV modules; finally, section five comments 
the conclusions of the work. 


2. Brief technology overview 


The main elements that constitute a HCPV generator are the 
solar cells, the receivers, the optical devices and the tracking 
system. A collection of solar receivers with electrically intercon- 
nected cells and with integrated optical devices is known as a 
HCPV module. Next, the main characteristics of these elements are 
described. 


2.1. Solar cells 


High concentration solar cells have specific characteristics 
different to those of conventional flat-plate solar cells. The HCPV 
cells are designed to extract more energy under a larger solar 
radiation. These cells have to be exclusively made of high quality 
semiconductor materials because of the significance of the con- 
version efficiency, a factor which is directly related to the quality 
of the crystal used, in which impurities must be avoided. 

An important characteristic of these solar cells is their 
enhanced thermal behavior. Under high temperatures, both Si 
and GaAs cells suffer a significant decrease of performance mainly 
because of the voltage reduction produced by a temperature 
increase. Nonetheless, this effect is less accurate at high concen- 
tration levels. 

Despite the fact that many different types of cells are being 
tested at laboratory level and although the reported efficiency has 
increased over 40% for multijunction cells [10], many of the 
concentration systems installed nowadays use silicon solar cells, 
providing efficiencies which hardly exceed 25% under concentra- 
tion. However, most of the current developments are based on 
triple-junction cells. These cells can be described as a stack of cells 
of different compositions with a progressively decreasing band- 
gap, so that the cell on the top of the stack absorbs higher energy 
photons than the cell in the middle, and this absorbs higher 
energy photons than the cell in the bottom. 


2.2. Receivers 


Solar cells have to be mounted on a receiver to extract the 
generated current and to dissipate and remove the heat resulting 
from sun power. Current extraction is made by thick wires or 
ribbons that are soldered to the cells, whereas heat removal is 
carried out by a heat sink added to the receiver. Materials like 
copper or aluminum are currently used for heat evacuation. Also, 
the solar receivers often include a bypass diode connected in 
parallel with the cell in order to avoid overheating of the cell in 
non-usual operating conditions such as shading. 


2.3. Optical devices 


The aim of the optical devices used in this technology is to 
concentrate radiation and to increase the luminous flux on the 
cells. Optical devices may be simple or consisting of primary and 
secondary optical element (SOE). Primary optical devices collect 
direct sun rays. Secondary optical devices receive the light from 
the primary ones and then change its properties by means of 
spectral filtering, light homogenization or light direction changing. 
In this respect, different optical devices are currently used: circular 
parabolic dishes, parabolic dishes with secondary, square flat 
Fresnel lenses, square flat Fresnel lenses with secondary, linear 
flat lenses, linear arched lenses or linear parabolic reflectors. 

Most of the systems currently designed are based on Fresnel 
lenses. This kind of lens reduces the amount of material required to 
concentrate the light by splitting the lens into a set of concentric 
annular sections known as Fresnel zones. The use of these zones 
allows keeping the required curvature without increasing the 
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thickness by means of adding discontinuities between them. Fig. 1 
shows an example of a solar receiver integrated with point-focus 
Fresnel lens. 


2.4. HCPV module 


A HCPV module is the smallest, complete, environmentally 
protected assembly of receivers and optical devices that is able to 
transform an input of unconcentrated solar radiation. Consequently, 
a typical HCPV module is made up of a group of cells, primary optics, 
secondary optics (optional) and housing components such as 


Primary { 
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Fresnel lens 


SOE 


Solar 
receiver 


Diode Solar cell Heat sink 


Fig. 1. Schematic of a HCPV solar receiver integrated with point-focus Fresnel lens. 
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Fig. 2. Schematic of a HCPV module based on point-focus Fresnel lenses. 
Source: IEC 62108 [2]. 
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interconnection and mounting. Fig. 2 shows an example of a typical 
HCPV module and its components. 


2.5. Tracking system 


Point-focus based HCPV modules must be always mounted on 
two-axis solar trackers, i.e. the modules must be always pointing 
the solar rays in order for the lenses to be able to focus the 
radiation on the small solar cell area. This kind of tracking system 
is more complex than linear tracking or static systems from a 
mechanical point of view but it presents the advantage of max- 
imizing the solar radiation capture along the daily trajectory. 


3. HCPV cells models 


In this section, five existing models that reproduce the behavior 
of a multijunction solar cell operating under high concentration 
are summarized. The presented models assume certain simplifica- 
tions in order to be applicable in engineering studies. Other more 
sophisticated models that divide the cell into hundreds or thou- 
sands of elementary cells are not included in this review [11,12]. 
These distributed models are difficult to implement and require 
high computational effort, so that they are not commonly used at 
the engineering level. 


3.1. Models based on the equivalent circuit of the multijunction cell 


A triple junction solar cell can be considered as composed of three 
series-connected p-n junctions. Each junction can be represented by 
an equivalent circuit model, being the most used models the Single 
Exponential Model (SEM) and the Double Exponential Model (DEM). 
Both models differ in the number of diodes that characterize the 
saturation current. By connecting the three subcells circuits in series, 
the equivalent circuit for the triple junction cell is obtained, Fig. 3. 

Taking into account the DEM model, the cell I-V curve can be 
mathematically expressed by the next set of equations 


V= Ù Vi 0) 


I= Ipp i—lo1itexp[(V; + IRsi)/m1iVr]-1} 
—lo2;{exp[(V; + IRs,)/M2;V7r] 
—1}-(V; + IRsi)/Rpi , i=1, 2, 3 (2) 


where the subscript i indicates the subcell (1=top, 2= 
middle, 3=bottom). Each subcell is characterized by a set of 


Rs,top | Rs,top | 
+ + ANV > 2 ~- s ry AM —— 2 
+ + + + 
Iph,top() Dtop X Rp,top = Vtop Iph,top() Dtop1 ¥ Dtop2 ¥ Rp,top = Vtop 
Rs,mid | Rs,mid | 
. . AW + + aT WV 
+ + 
Iph,mid) Dmid ¥  Rp,mid = Vmid V Iph,mid 8) Dmidt ¥  Dmid2 ¥ Rp,mid = uma Y 
Rs,bot | Rs,bot | 
s * AW Di + aT AWN 
= + 
iphibotÑ) Dbot ¥ Rp,bot $ Vbot Iphbot(Ñ) bot! ¥ Dbot2 ¥ Rp,bot $ Vbot 


Fig. 3. Equivalent circuit models for triple junction solar cells by using the Single Exponential Model (left) and the Double Exponential Model (right). 
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7 parameters: Ipp; (the photocurrent), Io1, and Io2, (the saturation 

currents), m4; and m32; (the diode ideality factors), Rs; (the series 

resistance) and Rp; (the shunt resistance). Vr is the semiconductor 

thermal voltage, Vr=kT/q, being k the Boltzmann constant, T the 

absolute temperature of the junction and q the electron charge. 
The equations for the SEM scheme are as follows: 


V= > Vi 8) 


T= Ipyi—loi{exp[(Vi + IRs i)/MiVr]—1} 
—(V; + IRsi)/Rpi , 1=1, 2, 3 (4) 


where each subcell is now characterized by 5 parameters: Ipp, (the 
photocurrent), Jo, (the saturation current), m; (the diode ideality 
factor), Rs; (the series resistance) and Rp; (the shunt resistance). 

In general, the model parameters depend on the operating 
conditions of the multijunction cell, such as concentration and cell 
temperature. However, as multijunction cells are often monolithic 
devices which incorporate only two electrical terminals, the sepa- 
rate measurement of the individual subcells is not possible. There- 
fore, the determination of the subcells parameters is very difficult. 

Some authors have proposed methods for extracting the sub- 
cells parameters, or at least some of them, by applying different 
assumptions and simplifications. We can find several extraction 
methods based on both the SEM scheme [13-18] and the DEM 
scheme [17,19-23]. However, it is remarkable that these studies 
often estimate parameters only at particular operating conditions 
and not for the whole range of possible conditions, so that the 
application of these models is sometimes limited. 


3.2. Model of G.S. Kinsey et al. (2008) 


The model proposed by G.S. Kinsey et al. [24,25] allows 
calculating the maximum power (Pypp) of a triple junction solar 
cell from its operating conditions: concentration ratio, X, cell 
temperature, Tc, and incident spectrum, F(A). First, three electrical 
parameters of the solar cell as a whole are calculated: the open- 
circuit voltage, Voc, the fill factor, FF, and the short-circuit current 
density, Jsc. Once these parameters have been determined, the 
Pmpp is obtained by: 


Pumpr = FEJscAcet Voc (5) 


Ace being the solar cell active area. 

The Voc and FF are assumed to be independent of the incident 
spectrum. While this assumption is adequate for the Voc [26,27], 
the FF variations with spectrum can have a significant impact on 
performance [28,29]. However, in this model, these variations are 
neglected. The Voc is obtained from the reference Voc at one sun 
and 25 °C of cell temperature by applying a logarithmic correction 
with concentration and a linear correction with temperature. 
On the other hand, the FF is obtained from the reference FF by 
only applying a linear correction with temperature. 

The model requires knowing the External Spectral Response of 
the component subcells (ESR;) as a function of the wavelength, 4, 
at different cell temperatures. These functions must be measured 
with the help of a solar simulator. A solar simulator usually 
consists of three major components: (1) light source(s) and 
associated power supply; (2) any optics and filters required to 
modify the output beam and (3) the necessary controls to operate 
the simulator, adjust irradiance, etc. [30]. 

The model equations for calculating the short-circuit current 
density are: 


eee | ESRtop(A. Tc)E(A)AA © 


JPHmid = | ESR mial, Tc )E()GA (7) 


Jsc = Min pytop-JpHmia} (8) 


where Jpytop, JpHmia are the photocurrents of the top and the 
middle subcells. Eq. (8) considers that the short-circuit current of 
the cell is limited by the smaller photocurrent of the top and 
middle subcells, i.e. the bottom subcell never limits the current. 

The model requires knowing several parameters that charac- 
terize the multijunction cell, being the most important the ESR of 
the top and the middle subcells. These values are obtained with 
the help of the HIPSS solar simulator from Spectrolab [31]. In 
addition, the use of a spectroradiometer is required to measure the 
incident spectrum. Alternatively, the spectrum can be estimated 
from atmospheric parameters by using simulation software such 
as SMARTS [32,33]. 


3.3. Model of C. Dominguez et al. (2010) 


The direct measurement of the spectrum is difficult and 
requires the use of a spectroradiometer, which is not always 
available. Another option is to use isotype cells. These sensors 
have the same spectral response than the component subcells of 
the multijunction device. This way, the photocurrents measured by 
the isotype cells under the same spectrum than the multijunction 
cell allow knowing the photocurrents of the component subcells. 

The model proposed by C. Dominguez et al. [16] is based on the 
use of isotype cells. It is a translation model, which calculates the 
points of the I-V curve at the desired operating condition from the 
points of the I-V curve at the reference condition. Three transla- 
tions are applied to each point of the reference I-V curve: (1) a 
translation of the current coordinate to the desired value of 
concentration; (2) a translation of the voltage coordinate to the 
desired values of concentration and incident spectrum; (3) a 
correction of the voltage coordinate to take into account the 
change in the cell operating temperature. This way, the model 
calculates the I-V characteristics as a function of irradiance, 
spectrum and cell temperature. 

Translation equations are analytically obtained by considering 
the SEM equivalent circuit model for multijunction cells (Fig. 3). 
In this scheme, the three shunt conductances are neglected and 
the three series resistances are lumped into a single effective 
series resistance, Rs. For instance, the translation of the voltage 
coordinate due to the change in concentration and spectrum is 
represented by the following equation 


k(Tcref + 273.15 
V(Tcret) = Vref A cret q ) 


id 
{mumo Hn (€E ot, ame ct) Tep jp mu 
PHtop,ref ‘ref PHmid,ref ‘ref 


= mbot 
(Fitouts ter) | jag DR (9) 


Ipubot,ref Tet 


(Vief, Iref) being the considered point of the reference I-V curve, 
V(Tcref) the translated voltage, Ipyirer the photocurrents of the 
component subcells at reference conditions, X the desired con- 
centration, m; the diode ideality factors of the component subcells 
and myy an effective ideality factor for the multijunction device. 
The factors F; are defined for each subcell as 


Ipi /Iphi ref 
—— - (10) 
min ({IpHtop. [PHmid> Tprbot}/MiN{Ipytop,ref» Tpumid,ref> /pHbot ret } 


Fi 


i.e. they account for the photocurrents of the subcells, Ipp;, as 
measured by the isotype cells. 
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The model depends on some parameters that characterize the 
multijunction device. Some of them can be found in semiconductor 
textbooks but some other must be fitted by using simulation tools 
that minimize the error between experimental data and modeled 
data. A solar simulator is required to obtain the experimental data. 


3.4. Model of E.F. Fernández et al. (2013) 


The aim of the model proposed by E.F. Fernandez et al. [34,35] 
is to obtain a model based on easily adjustable parameters. 

The measurements of several multijunction cells showed that 
the Ge bottom subcell never limits the short-circuit current of the 
multijunction device and has a small influence on the rest of the 
electrical parameters under normal operating conditions. There- 
fore, the behavior of the cell can be modeled by considering 
only two “equivalent subcells” that verify the equation of the 
SEM model. 

The short-circuit currents of the two equivalent subcells (Isc1, 
Isc) are supposed to be equal than the photocurrents of the top 
and the middle real subcells. These photocurrents must be 
measured with isotype cells or calculated from the incident 
spectrum as measured by a spectroradiometer if the ESR of the 
subcells is known. 

On the other hand, the open-circuit voltages of the equivalent 
cells (Voc1, Voc2) verify the condition that their sum is equal to the 
open-circuit voltage of the real cell. A logarithmic correction on 
concentration is used to calculate these open-circuit voltages from 
the Voc values at one sun. The objective of calculating these open- 
circuit voltages is to obtain the saturation currents of the equiva- 
lent cells (Io, Io2), which are computed by using well-known 
simplified expressions derived from the SEM solar cell model. 

The equation of the I-V curve of the multijunction cell is 
modeled by means of a series association of the two equivalent 
subcells, i.e. the voltages of these subcells are summed in order to 
obtain: 


V= mk(Tc + 273.15) n( “5 =) 41n a] -IRs (11) 
q To To2 


Rs and m being effective values for the series resistance and the 
diode ideality factor of the multijunction device. Finally, this LV curve 
is corrected taking into account an effective shunt resistance Rp. 

Because of the simplicity of the model equations, the imple- 
mentation of the numerical algorithm for obtaining the model 
parameters is easy. However, like in the models commented above, 
a set of measurements carried out in a solar simulator is required 
for adjusting these parameters. 


3.5. Model of Sandia National Laboratories (2004) 


The models commented above cannot always be applied in 
practice, because of the need of a solar simulator, the need of 
instruments such as spectroradiometer or isotype cells or the 
difficulty of implementing the algorithms for adjusting the model 
parameters. So, in some cases, the use of a simpler model could be 
justified in spite of a little loss in accuracy. The model developed at 
Sandia National Laboratories [36] was mainly conceived for 
characterizing modules and arrays of different technologies, but 
it also can be applied to high concentration solar cells. It does not 
require sophisticated instruments and the model parameters can 
be obtained from outdoor measurements. 

As it has been shown, the main problem for characterizing 
multijunction cells is that they are very sensitive to the incident 
spectrum. In general, the sunlight spectrum varies as a function of 
air mass, clouds, turbidity due to aerosol effects and precipitable 
water vapor [28]. However, the main contribution to the spectrum 


change is air mass [37,38]. In the model of Sandia National 
Laboratories, the measurement of the sunlight spectrum is not 
required and the spectral influences are only quantified by the air 
mass parameter, as defined in [39]. 

The model equations are formulated as a function of the 
“effective irradiance” (Ber). In the case of high concentration cells, 
this parameter is defined as: 


Bef = [Gef (AM)]/G3 (12) 


Gg being the incident direct irradiance and Gg* the reference 
direct irradiance. f;(AM) is the spectral correction factor as defined 
in the standard ASTM E 973 [40], which is approximated by a 
fourth order polynomial function of the air mass. 

Five points of the cell I-V curve are calculated by the model: the 
three classical points (short-circuit, open-circuit and maximum 
power points) and two additional points: (Vx, Ix), the point at a 
voltage Vy=Voc/2, and (Vxx,/xx), the point at a voltage Vxx= 
(Vupp+Voc)/2. Fig. 4 indicates the position of these points on the 
cell I-V curve. As an example of the model equations, the 
expressions that allow calculating the Isc and Voc parameters are 
the following 


Isc = Berle + aisc(Tc-TQ)] a3) 


Voc = Vic + SINBet) + Byoc(Tc—T2) oa 


where the superscript “+” indicates parameters at reference con- 
ditions. In these equations, the influence of cell temperature is 
modeled with the temperature coefficients ajs- and Pyoc, where 
PBvoc is allowed to vary linearly with the effective irradiance. 
Coefficient 5 is the product of the effective diode ideality factor, 
m, and the thermal voltage Vz. 

Every model parameter can be obtained by regression analysis 
from outdoor monitorized data, what is an advantage with respect 
to other models. The main disadvantage of this model is that it 
does not take into account ambient factors that affect the sunlight 
spectrum, such as clouds, aerosols and precipitable water vapor. 
The model is oriented to clear days with clean atmosphere. 
However, it is remarkable that these conditions are prevalent at 
the locations where HCPV systems are installed. 

Table 1 summarizes the main characteristics of the referenced 
models for high concentration photovoltaic cells, including the 
required inputs and advantages of each one. 


(0,Isc) (Vx, Lx) 
(Vmpp,Impp) 
i= 
g 
5 
O 
0 
0 Vx=Voc/2 Vmpp Vxx= Voc 
Voltage (Vmpp+Voc)/2 


Fig. 4. Points of the I-V curve calculated by the model of Sandia National 
Laboratories. 
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4. Models for estimating the maximum power of a HCPV 
module 


This section summarizes the state-of-the-art of models for 
HCPV modules characterization. These models take into account 
different atmospheric parameters that affect the behavior of the 
module, such as direct normal irradiance (DNI), air mass (AM), air 
temperature (Tir), wind speed (WS), etc., and some of the models 
use the mean cell temperature of the module (Tc) as input. 


4.1. Model of the standard ASTM E 2527 


The standard ASTM E 2527 [41] defines a very simple way of 
estimating the Pypp of a high concentration module. It uses as 
inputs the direct normal irradiance (DNI), the air temperature (Tair) 
and the wind speed (WS). It is based on a single equation: 


Pumpr = DNI(c1 + C2DNI 4 C3T air 4 c4WS) 


(15) 


where c1, C2, C3 and c4 coefficients are obtained by regression 
analysis of outdoor monitorized data. As can be seen, the model 
does not take into account the spectral effects that significantly 
affect the behavior of modules containing multijunction solar cells. 


4.2. Model of G. Peharz et al. (2011) 


The model of G. Peharz et al. [42] does take into account the 
spectral influences. It is based on the following equation: 


Pupp = CpntDNI + Cz2Z° + C7Z + CreTc + Coftset (16) 


where the Z parameter represents the mismatch between the 
photocurrents of the top and middle subcells with respect to the 
reference AM1.5d spectrum. The value of Z is obtained from the 
photocurrents as measured by isotype cells [43]. A value of Z>0 
corresponds to a blue spectrum, where the limiting subcell is the 
middle subcell, while a value of Z<O corresponds to a red 
spectrum, where the limiting subcell is the top subcell. In order 
to quantify the value of Z, some parameters of the solar cell must 
be measured in a solar simulator. Cpn, Cz, Cz, Cre and Coftset 
coefficients are obtained by regression analysis from outdoor 
monitorized data. 


4.3. Model of M. Steiner et al. (2012) 


The model proposed by M. Steiner et al. [44] is a complete 
model that combines the spectrum simulation software SMARTS 
[32,33] with ray tracing software for modeling the optics of the 
Fresnel lenses [45] and a SPICE network model for modeling the 
I-V characteristics [46,47]. The model requires a complete set of 
instruments in order to be applied because it takes into account 


Table 1 


many factors that affect the module behavior. Beyond a pyrheli- 
ometer and weather station equipment, it requires isotype cells in 
order to measure the Z parameter, multifilter rotating shadowband 
radiometer for the determination of Aerosol Optical Depth (AOD) 
and Precipitable Water (PW) and tracking accuracy sensor for 
measuring the alignment of the tracker to the sun. The model 
allows calculating the Pyypp of high concentration modules with a 
high level of accuracy, being its main disadvantages the difficulty 
of implementing and the equipment requirements. 


4.4. Model of A.J. Rivera et al. (2013) 


Due to the fact that the relation between atmospheric para- 
meters and module output maximum power is complex and non- 
linear, the use of artificial neural network-based models has been 
proposed. This approach has been followed by AJ. Rivera et al. [48]. 
In particular, a cooperative-competitive hybrid algorithm for radial 
basis function networks was implemented. 

The model uses the Average Photon Energy (APE) [49,50] in 
order to quantify the spectral influences on HCPV module max- 
imum power. This is a single value in eV that characterizes the 
shape of the spectrum, so that the higher APE, the higher the blue 
content of light. A spectroradiometer is required for obtaining 
this input. 

The model takes into account the direct normal irradiance 
(DNI), the air temperature (Tair), the wind speed (WS) and the APE. 
Coefficients of the neural network are obtained from outdoor 
monitorized data. 


4.5. Model of F. Almonacid et al. (2013) 


The solar research group of Jaén University has wide experi- 
ence on the application of artificial neural networks in the 
photovoltaic field [51-56]. Taking into account this experience 
and from a similar point of view than the previously commented 
model, the model proposed by F. Almonacid et al. [57] tries to 
characterize the relation between atmospheric parameters and 
module output maximum power through artificial neural net- 
works. In this case, a feed-forward neural network trained with 
the Levenberg-Marquardt back-propagation algorithm was used. 

The model takes into account the spectral influences through 
easily measurable parameters: air mass (AM) and precipitable 
water (PW). This way, it avoids the need of spectroradiometer. 
Inputs of the model are the direct normal irradiance (DNI), the air 
temperature (Tair), the wind speed (WS), the AM and the PW. 
Coefficients of the neural network are obtained from outdoor 
monitorized data. 
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Models: Equivalent circuit Kinsey Dominguez Fernandez Sandia 
Required inputs 

Concentration ratio y y y y y 
Cell temperature y y y y y 
Incident spectrum y y 

Subcells photocurrents y y 

Air mass y 
Advantages 

Provides the whole I-V curve y y y 

Spectral effects quantified with high accuracy y y y y 

Explicit equations y y y y 
Algorithms for extracting the model parameters are easy to implement y y y 
It does not require knowing in detail the incident spectrum y y y 
It does not require indoor characterization of the solar cell y 
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4.6. Model of E.F. Fernández et al. (2012) 


The model of E.F. Fernandez et al. [35] uses air mass (AM) in 
order to quantify the spectral influences on HCPV module max- 
imum power. Therefore, it avoids the need of spectroradiometer or 
isotype cells. The aim is to have a model based on easily measur- 
able parameters. The main equations of the model are as follows: 


Pypp = A; DNI +AT air for AM <2 (17) 


Pypp = A1 DNI + AdT air + A3AM for AM > 2 (18) 


These equations consider that, for large values of AM, Pypp is 
mainly influenced by DNI, Tay and AM, while for small values of 
AM, the influence of AM can be neglected [58,59]. Coefficients Aj, 
Az and A, are obtained by regression analysis from outdoor 
monitorized data. At present, the model has been only evaluated 
at the south of Spain. 


4.7. Model of Sandia National Laboratories (2004) 


The models above only calculate the maximum power of the 
module. When more information about the device I-V curve is 
required, an alternative is to use the model of Sandia National 
Laboratories [36,60-63], which calculates five points of the I-V 
curve as shown in Fig. 4. 

Equations of the model are formulated as a function of the 
“effective irradiance”, which, for the case of a HCPV module, is 
defined as it is shown in Eq. (12). 

The model equations that allow calculating the maximum 
power of a HCPV module are: 


fi (AM) = dp + d1 AM + a)AM? + 3AM? + a4AMÎ (19) 
Ber = [DNI f(AM)]/DNI* (20) 
6=[mk(Tc + 273.15)]/q (21) 
Pympp = Pvmpp o + Mvmpp(1—Ber) (22) 
Impp = (CoBer + C1Bep)l pp + Aimpp(Tc-Te)] (23) 


Vmpr = Viupp + C2Nsdln(Ber) + C3Ns[5In(Ber)]? + Bympp(Tc-Te) 24 


Table 2 


Pumpr = Impp Vee (25) 


Eq. (19) approximates the spectral correction factor, fi(AM), 
as defined in the standard ASTM E 973 [40], by a fourth order 
polynomial. Eq. (20) calculates the effective irradiance (Ber), that is, 
the irradiance at which the cells actually respond. This irradiance 
is the direct normal irradiance, DNI, corrected with the spectral 
correction factor, f;(AM), and normalized to the reference irradi- 
ance, DNI*. Eq. (21) defines the 6 parameter, which is the product 
of the effective ideality factor of the multi-junction cell (m) and 
the thermal voltage. The thermal voltage is obtained from the 
Boltzmann constant (k), the electron charge (q) and the cell 
temperature (Tc). Eq. (22) allows the determination of the tem- 
perature coefficient Sympp. This coefficient is used afterwards for 
quantifying the effect of temperature on module maximum power 
point voltage. The coefficient is expressed as a linear function of 
the effective irradiance, i.e. it is allowed to vary with the concen- 
tration ratio. Eqs. (23) and (24) calculate the module maximum 
power point current (Impp) and voltage (Vmpp) from their values at 
reference conditions (Iypp*, Vmpp*). The reference conditions are 
defined in the model by DNI*, Tc* and AM1.5. aimpp is the 
temperature coefficient for Impp and Ns is the number of cells in 
series for the module. Finally, Eq. (25) obtains the maximum 
power of the HCPV module. 

As it has been shown, some of the referenced module models 
use the mean cell temperature of the module (Tc) as input. The 
direct measurement of this parameter is very difficult, due to the 
special features of these kinds of devices. Because of this, indirect 
methods for estimating Tc have been proposed [42,6469]. 

Table 2 summarizes the main characteristics of the referenced 
models for HCPV modules, including the required inputs and 
advantages of each one. 

Table 3 presents a quantitative comparison of the referenced 
models for HCPV modules characterization. Four models have 
been evaluated at the Center of Advances Studies in Energy and 
Environment of Jaén University (ASTM E 2527, Almonacid, 
Fernandez and Sandia). The rest of the models have been eval- 
uated by their respective authors and results can be found in the 
papers that describe the models. With respect to the values shown 
in Table 3, it is important to remark that a direct comparison of the 
provided values is not possible because the different authors use 
different measurement and evaluation procedures. Furthermore, 
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Models: 


Required inputs 
Direct normal irradiance 
Ambient temperature 
Wind speed 
Precipitable water 

Cell temperature 
Incident spectrum 
Subcells photocurrents 
Average photon energy 
Aerosol optical depth 
Air mass 

Tracking error 


LAS 


Advantages 

It quantifies spectral effects 

Spectral effects quantified with high accuracy 

It does not require the cell temperature 

It does not require spectroradiometer 

It does not require indoor characterization 

It provides electrical parameters other than Pmpp 

Easy to implement 

It does not require advanced knowledge on artificial neural networks 
Low computational effort 
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Table 3 
Errors between measured and calculated maximum power of the models for HCPV 
modules as provided by different authors. 


Model Parameter Value 

ASTM E 2527 RMSE 4.6%* 

Peharz Absolute RMSE 1.3, 1.2, 1.6 and 0.6 W** 
Steiner NRMSE Between 2.5 and 3.2%*** 
Rivera MSE 65° 

Almonacid RMSE 2.1%* 

Fernandez RMSE 3.2%* 

Sandia RMSE 3.4%" 


* Mean value of the root mean square error for two HCPV modules measured 
during 2011 and 2012 at the Center of Advanced Studies in Energy and Environ- 
ment of Jaén University (south of Spain). 

** Values of the absolute root mean square error for four HCPV modules 
measured in different time periods from 2009 to 2012 at the Fraunhofer ISE 
(Freiburg, Germany) [42]. The nominal power of the modules at 850 W/m? of direct 
normal irradiance, 25 °C of cell temperature under the incident spectrum AM 1.5d 
ASTM G173-03 is 54.0, 50.1, 44.0 and 15.7 W respectively. 

** Range of the obtained normalized root mean square error for five HCPV 
modules measured from November 2011 to March 2012 at the Fraunhofer ISE 
(Freiburg, Germany) [44]. The errors are normalized to the nominal power of each 
module. 

** Value of the mean square error for one HCPV module measured from 
November 2011 to March 2012 at Jaén University (south of Spain) [48]. 


the reported errors are expressed by different parameters in each 
case. However, the table can be illustrative of the order of 
magnitude of the errors of the models. 


5. Conclusions 


Nowadays, we can find in the scientific literature several models 
for the electrical characterization of HCPV cells and modules. These 
models are required for the design, monitoring and energy prediction 
of HCPV systems. A review of the existing models has been presented 
with the aim of helping the photovoltaic professionals and research- 
ers in the choice of the most suitable model for each specific 
application. 

Multijunction solar cell models calculate the relevant electrical 
parameters of the cell as a function of irradiance, spectrum and cell 
temperature. In order to evaluate the spectral influences, spectro- 
radiometer, isotype cells or a solar simulator can be required. 
However, we can find other simplified models that require less 
sophisticated instruments in spite of a little loss of accuracy. 
Numerical algorithms are often needed in order to adjust the model 
parameters. These algorithms are not always easy to implement. 

HCPV modules models calculate the relevant electrical para- 
meters as a function of ambient parameters such as direct normal 
irradiance, air temperature or wind speed. The spectral influences 
can be characterized in different ways (Z, APE, AM) and different 
equipments are required to obtain these spectral parameters. 
Every model requires outdoor measurements in order to obtain 
the model coefficients by means of regression analysis. Some of 
these models are easy to implement and do not require very 
sophisticated instruments to be applied so that they are good 
candidates to be used in the industry, in particular the ASTM E 
2527 [41], the model developed by Fernandez et al. [35] and the 
model developed by Sandia National Laboratories [36]. 
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